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Abstract 

At the Large Hadron Collider, we prove the feasibility to detect pair production of the 
lightest CP-even Higgs boson h of Type II 2-Higgs Doublet Models through qq( ) — > 
qqC^hh (vector-boson fusion). We also show that, through the hh — > 46 decay channel in 
presence of heavy-flavour tagging, further exploiting forward/backward jet sampling, one 
has direct access to the Xuhh triple Higgs coupling - which constrains the form of the 
Higgs potential. 
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1 Introduction 



If only a light Higgs boson (with mass Mh ^140 GeV) is found at the Large Hadron Collider 
(LHC), it may be difficult to tell whether it belongs to the Standard Model (SM) or indeed 
a model with an enlarged Higgs sector. For example, in the case of a CP-conserving Type 
II 2-Higgs Doublet Model (2HDM) [lj-[4:jf[ possibly in presence of minimal Supersymmetry 
(SUSY) - the combination of the two yielding the so-called Minimal Supersymmetric Standard 
Model (MSSM) - this happens in the so-called 'decoupling region', when M H , Ma, M H ± 3> M h , 
for suitable choices of the other MSSM and 2HDM parameters, where - for the same mass - the 
h couplings to ordinary matter in the SM are the same as in both the 2HDM and MSSM. Even 
in these conditions, however, it has been proved that one could possibly establish the presence 
of an extended Higgs sector by determining the size of the trilinear Higgs self-coupling X hhh [5]. 

If the extended model is not in a decoupling condition, then it is generally possible to 
establish the presence of additional Higgs signals, H, A and/or [SI El- However, even when 
this is the case, it may be difficult to distinguish, e.g., between a generic Type II 2HDM and 
the MSSM (unless, of course, one also detects the SUSY partners of ordinary matter and Higgs 
bosons). In fact, despite there exist well establish spectra among the four different masses in the 
MSSM (for fixed, say, Mh and tan /3, the ratio of the vacuum expectation values of the two Higgs 
doublets in either model), it may well be possible that the additional 2HDM parameters arrange 
themselves to produce an identical mass pattern. However, such a degeneracy between the two 
models would not typically persist if one were able to also measure certain Higgs couplings, 
chiefly those among the Higgs bosons themselves (involving two or more such particles). In 
fact, while the measurement of only two among the four Higgs boson masses (Mh, Mh, Ma and 
Mh± ) - or, alternatively, one such masses and tan (3 - would fix (at tree-level) all Higgs masses 
and couplings in the MSSM, this is no longer true in a generic Type II 2HDM [1], because of the 
freedom in selecting the free additional parameters. For example, the general CP-conserving 
Type II 2HDM that we are going to consider can be specified uniquely by seven parameters: 
Mh, Mh, Ma, M h ±, (3, a (the mixing angle between the two CP-even neutral Higgs states) and 
A 5 (see eq. ([2]) later on). It may then happen that the first six of these are measured and found 
to agree with the MSSM pattern, but one would still need to measure A5 to verify that it is the 
Higgs sector of the MSSM that is present. One way to do so would be by measuring trilinear 
Higgs self-couplings, such as Xhhh and \hmi- Alternatively, the measurement of the latter two 
couplings would constitute a test of the MSSM relations if one knew Mh and tan /3 but not a. 

In this paper, we make the assumption that only one parameter is known, Mh, as may 
well happen at the LHC after only a h resonance is detected. We further imply that all 

2 Of the initial eight degrees of freedom pertaining to the two complex Higgs doublets, only five survive as 
real particles upon Electro- Weak Symmetry Breaking (EWSB), labelled as h, H, A (the first two are CP-even or 
'scalars' (with Mj, < Mh) whereas the third is CP-odd or 'pseudoscalar') and i? ± , as three degrees of freedom 
are absorbed into the definition of the longitudinal polarisation for the gauge bosons Z and W , upon their 
mass generation after EWSB. 



(potential) SUSY states are much heavier than ordinary particles (with the possible exception 
of the lightest SUSY particle, see footnote 4), thus effectively a decoupled MSSM setup. Under 
these circumstances, we then ask ourselves the following question. While trying to establish the 
presence of additional (single) heavy Higgs signals, which would then unmistakably distinguish 
between the SM and a scenario with an extended Higgs sector, would it also be possible to gather 
information on Higgs self-couplings from signatures involving two light Higgs bosons, hence by 
studying channels involving h pair production, thereby possibly also distinguish between, e.g., 
a generic Type II 2HDM and the MSSM? 

It is the purpose of this paper to show that this is the case, so long that enough luminosity 
can be accumulated at the LHC, also in view of the Super-LHC (SLHC) option [Sj. We will 
illustrate how we have come to this conclusion, i.e., after investigating the process [H] 

qq^ — > qq^ ^hh (vector — boson fusion), (1) 

with referring to any possible (anti)quark flavour combinations^. The relevant Feynman 
diagrams corresponding to process ([1]) in both the MSSM and 2HDM considered here can be 
found in Fig. [TJ In our selection analysis, we will resort to the extraction of two h — > bb 
resonances, in presence of the following signature: 

• 'four 6-quark jets and two forward/backward-jets'. 

This signature was already considered in Ref. [5] in the SM context (from which we will import 
some of the results). 

Our paper is organised as follows. In the next section, we outline the computational proce- 
dure. Sect. 3 presents our numerical results and discusses these in various subsections. Sect. 4 
contains our conclusions. 

2 Calculation 

We have assumed y/s = 14 TeV for the LHC energy throughout. Our numerical results are 
obtained by setting the renormalisation and factorisation scales to 2M/, for the signal while for 
the QCD background we have used the average jet transverse momentum (pj, = YliPTj/ n )- 
Both Higgs processes and noise were estimated by using the Parton Distribution Function 
(PDF) set MRST99(COR01) [IS) . While the background calculations were based on exact tree- 
level Matrix Elements (MEs) using the ALPGEN program [19] . all signal rates were obtained 

3 The gluon-gluon production mode [10] was considered in Refs. [11] and [12] (see also [13]), and later on 
[141115]. where - despite significant kinematic differences exist between signal and QCD noise - it was eventually 
shown that the extraction of the gg — > hh — > bbbb signal is essentially impossible at the (S)LHC because of the 
overwhelming QCD noise, both reducible and irreducible. Recently, encouraging results on the cross-section for 
multi-Higgs boson production in the gluon-gluon production mode has been obtained in models beyond the SM 
and MSSM [16]. The possibility of using Higgs boson pair production more generally to access trilinear Higgs 
couplings has also been studied on the level of total cross-sections in [T7] . 



through the same level of accuracy via programs based on the HELAS subroutines [20] - for 
the computation of the MEs - and VEGAS [21] or Metropolis [22] - for the multi-dimensional 
integrations over the phase space. As for numerical input values of SM parameters, we adopted 
the ALPGEN defaults. 

Concerning the MSSM setup, the two independent tree-level parameters that we adopt are 
Ma and tan/?. Through higher orders, we have considered the so called 'Maximal Mixing' 
scenario (X t — A t — fi/ tan/3 = v^^susy) [23], wherein we have chosen for the relevant SUSY 
input parameters: /i = 200 GeV, A b = 0, with M SU sy = 5 TeV, the latter - as already intimated 
- implying a sufficiently heavy scale for all sparticle masses, so that these are not accessible at 
the LHC and no significant interplay between the SUSY and Higgs sectors of the model can 



take plact^- Masses and couplings within the MSSM have been obtained by using the HDECAY 
program [2"3] . 

Before giving the details of the 2HDM setup we are using, let us recall the most general 
CP-conserving 2HDM scalar potential which is symmetric under $1(2) — > —$1(2) up to softly 
breaking dimension-2 terms (thereby allowing for loop-induced flavour changing neutral cur- 
rents) [I], 

V = m?!$t$i + m|,$£$ 2 - {m? 2 $}$ 2 + h.c.) + ($I$i) 2 + ^X 2 (H^f + 

+A 3 ($l$i) ($J$ 2 ) + A 4 ($1$ 2) + {^A 5 ($ \<S> 2 ) 2 + h.c.} . (2) 

In the following, the parameters m n , m 2 2, rn 12 , X±, A 2 , A 3 and A 4 are replaced by v, M h , M H , 
Ma-i Mh±, (3 and a (with v fixed). Hence, as intimated already, the CP-conserving 2HDM 
potential is parameterised by seven free parameters. Notice that from the scalar potential all 
the different Higgs couplings needed for our study can easily be obtained. (See [21 [3] for a 
complete compilation of couplings in a general CP-conserving 2HDM.) 

In our 2HDM, we will fix M% and Mu to values similar to the ones found in the MSSM 
scenario we are considering, by adopting three different setups: 

1. M h = 115 GeV, Mjj = 300 GeV, 

2. M h = 115 GeV, M H = 500 GeV, 

3. M h = 115 GeV, M H = 700 GeV. 

We always scan over the remaining parameters in the ranges 

— 7r/2 < a < it/2, 
-An < A 5 < 47r, 
< tan/3 < 50, 



4 The only possible exception in this mass hierarchy would be the Lightest Supersymmetric Particle (LSP), 
whose mass may well be smaller than the lightest Higgs mass values that we will be considering. However, we 
have verified that invisible h decays (including the one into two LSPs) have negligible decay rates. 



100 GeV < M A < 1000 GeV, 



100 GeV < M H ± < 1000 GeV. 



In order to accept a point from the scan we also check that the following conditions are 
fulfilled: the potential is bounded from below, the A, fulfill the tree-level unitarity constraints 
of [25] and yield a contribution to |Ap| < 10~ 3 . In short the unitarity constraints amounts 
to putting limits on the eigen values of the S matrices for scattering various combinations of 
Higgs and electroweak gauge bosons. We have followed the normal procedure pQ of requiring 
the J = partial waves (ao) of the different scattering processes to fulfill |Re(ao)| < 1/2, which 
corresponds to applying the condition that the eig Ay^. ± of the scattering matrices (or 

more precisely 16ttS) fulfill |Ay^ ± | < 8tt |26j . In other words we allow parameter space points 
all the way up to the tree-level unitarity constraint |Re(ao)| < 1/2. In order to investigate 
the sensitivity to this upper limit we will also report results as a function of the value of the 
maximal eigenvalue, A max . The spectrum of masses, couplings and decay rates in our 2HDM 
is the same as in Ref. |27j . obtained by using a modification of HDECAY [21] (consistent with 
a similar manipulation of the program used in Ref. [2B])- For each accepted point in the scan 
the partial decay rates for the different Higgs bosons are then calculated using HDECAY and 
also taking possible additional partial widths of the H into account. 

While the parameter dependence of the MSSM Higgs sector renders the computation of the 
tree-level MSSM cross-sections rather straightforward (as the latter depends on two parameters 
only, Ma and tan f3), the task becomes much more time-consuming in the context of the 2HDM. 
In order to calculate the cross-sections in this scenario, they are schematically written as a 
combination of couplings and kinematic factors in the following way: 



where all the explicit dependence on a, /3, \Hhh and Xhhh is contained in the couplings gc. 
gi = sin 2 (P - a), g 2 = cos 2 (/? - a), g 3 = cos(/3 - a)X Hhh , g 4 = sin(/? - a)\ hhh , and g 5 = 1, 
whereas the dependence on masses and other couplings is in the factors 



Note that the sum over subamplitudes M t also contains all interference terms and that colour 
factors etc. are included properly^ The <Xy are then calculated numerically for fixed masses. 
We can then get the cross- sect ion in an arbitrary parameter space point by multiplying the 
kinematic factors with the appropriate couplings. However, there is a slight complication since 
the kinematic factor for the H — > hh contribution depends on the width T# if there is a s- 
channel resonance and the width in turn depends on the couplings. In this case the kinematic 
factor scales as 1/Yh which is accounted for by assuming a fixed value for the width when the 

Here, Zi refers to the Zi symmetry, Y is the hypercharge, and a is the total weak isospin 
6 We have carefully verified the integrity of our procedure. 






(4) 



kinematic factor is calculated and then rescaling the result with the true width when calculating 
the contribution to the cross-section. Another complication is the dependence of the kinematic 
factors on the Higgs masses, Ma and Mff±. The contributions of main interest, which contain 
the Xhhh an d Xhhh couplings, only depend on these masses indirectly through the unitarity 
constraints. At the same time there are other contributions to the cross-section which depend 
explicitly on these masses. However, these contributions are very small in the parts of parameter 
space of interest and can thus be safely neglected. 

3 Results 

In our investigation of the emerging hadronic final state, we will assume that 6-quark jets 
are distinguishable from light-quark and gluon ones and neglect considering 6-jet charge de- 
termination. Finite calorimeter resolution has been emulated through a Gaussian smearing 
in transverse momentum, p T , with (cr(p T )/p T ) 2 = (0.60/y^) 2 + (0.04) 2 , for all jets. The 
corresponding missing transverse momentum, p™ lss , was reconstructed from the vector sum of 
the visible momenta after resolution smearing. Finally, in our parton level analysis, we have 
identified jets with the partons from which they originate and applied all cuts directly to the 
latter, since parton shower and hadronisation effects were not included in our study. 

3.1 Inclusive Signal Results 

In this section, after a preliminary analysis of the Higgs mass and coupling spectra in the 
MSSM and a general Type II 2HDM, we will start our numerical analysis by investigating the 
model parameter dependence of the Higgs pair production process in (TjQ) at fully inclusive level, 
in presence of the decay of the latter into two bb pairs, with the integration over the phase 
space being performed with no kinematical restrictions. This will be followed by an analysis of 
the production and decay process pertaining to the Higgs signal of interest at fully differential 
level, in presence of detector acceptance cuts and kinematical selection constraints. Finally, we 
will compare the yield of the signal to that of the corresponding background and perform a 
dedicated signal-to-background study including an optimisation of the cuts in order to enhance 
the overall significance. We will treat the MSSM and 2HDM in two separate subsections. 

3.1.1 MSSM 

As representative of the low and high tan j3 regime, we will use in the remainder the values of 3 
and 40. We have instead treated Ma as a continuous parameter, varying between 100 and 700 
GeV or scB Before proceeding with the numerical analysis of the signal, it is worthwhile to 
investigate both the Higgs mass and coupling dependence in the MSSM with respect to the two 
input parameters Ma and tan {3. This is done in Figs. [2] and El respectively. In the latter, we 

7 Values of Ma below 90 GeV or so are actually excluded by LEP for the lower tan/3 value: see [25] . 



study the case of (CP-even) MSSM Higgs boson couplings to gauge bosons (denoted by Ghvv 
and Ghvv) j wherein V refers to either a or a Z. In the same figure, the symbol </> refers to 
the SM Higgs boson, with mass identical to that of the lightest MSSM Higgs state (M^ = Mh). 
While the pattern of masses has been well established in past literature, it is interesting to notice 
here that the product of the MSSM couplings entering process (1) is always smaller than in 
the SM case. However, in the MSSM, resonance enhancements can occur (such as in H — > hh), 
so that the actual MSSM production rates can in some cases overcome the corresponding SM 
ones (for M <t> = M h ). 

Fig. H] presents the fully inclusive MSSM cross-section for the process of interest, as defined 
in (JT]), times (effectively) BK(hh — ► bbbb). The shape of the curves is mainly dictated by the 
interplay between phase space (see Fig. [2]) and coupling (see Fig. [3]) effects, with the exception 
of the region Ma ~ 220 GeV and tan/? = 3, where the onset of the H — > hh resonance is 
clearly visible. Cross-sections are generally sizable, particularly at low tan/3. The displayed 
rates however coincide to the ideal situation in which all final state jets are detected with unit 
efficiency and the detector coverage extend to their entire phase space, so that they only serve 
as a guidance in rating the phenomenological relevance of the process discussed. 

A more realistic analysis is in order, which we have performed as follows. The four 6-jets 
emerging from the decay of the hh pair are accepted according to the following criteria: 

p b T > 30 GeV, \r] b \ < 2.5, AR bb > 0.7, (5) 

in transverse momentum, pseudorapidity and cone separation, respectively. Their tagging 
efficiency is taken as e& = 50% for each b satisfying these requirements, e b = otherwise^. In 
addition, to enforce the reconstruction of the two Higgs bosons, we require all such 6's in the 
event to be tagged and that at least one out of the three possible double pairings of 6-jets 
satisfies the following mass preselection: 

(m bl , h - M h f + {m HM - M h f < 2 a 2 m , (6) 

where a m = 0.12 M h . We further exploit 'forward/backward-jet' tagging, by imposing that the 
non-6-jets satisfy the additional cuts 

^fwd/bwd > 2Q GeV; 2 5 < ^ fwd < 5j _ 2 5 > ^ bwd > _^ (7) 

Tab.dJshows the rates of the signal after the implementation of the constraints in eqs. ©-(JTj) 
(hereafter, referred to as 'acceptance and preselection cuts' or 'primary cuts'). While our process 
does yield non-negligible rates after the latter, it turns out that it is of no phenomenological 
relevance, even assuming very high luminosity. Firstly, in view of the fact that 6-tagging 
efficiencies are not taken into account in this table: for the '46-jet' tagging option, one should 
multiply the numbers in Tab. [T]by ef, that is, 1/16. (Alternative approaches requiring a lesser 



8 Here and in the remainder, the label b refers to jets that are 6-tagged while j to any jet (even those 
originating from 6-quarks) which is not. 
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Table 1: Cross-sections for Higgs pair production via vector-boson fusion, process (JTJ), after 
Higgs boson decays (relevant BRs are all included) and the acceptance and preselection cuts 
defined in (j5])-([7j), for two choices of tan/? and a selection of Ma values, assuming the MSSM in 
Maximal Mixing configuration (the corresponding values of Mh are also indicated in brackets). 
No 6-tagging efficiencies are included here. 



number of 6-jets to be tagged as such were not successful either.) Secondly, the background 
rates, after the same cuts in eqs. are always overwhelming the signal, despite our efforts 

in further optimising the cuts. For this reason, rather than dwelling upon the latter now, we 
postpone their discussion to the next subsection and simply conclude here that our channel is 
altogether inaccessible at both the LHC and SLHC in the context of the MSSM. 



3.1.2 2HDM 

As already alluded to, the parameter space of the general CP-conserving Type II 2HDM we 
are considering is quite large as it depends on seven unknown parameters. In order to get a 
feel for the dependence of the signal cross-section for the process qq^ — > qq^hh — > qq^bbbb 
we therefore present in Figs. [5] through [7] the results of our three selected scenarios, wherein we 
scan the allowed parameter space over 10000 randomly chosen points. (Note that similarly to 
the MSSM case we have included the BK(hh — > bbbb) but not any 46-jet tagging efficiency.) 

Comparing with the cross-sections in the MSSM the main differences are due to the follow- 
ing: 

• the triple Higgs couplinggj Xuhh and Xhhh are not related to the gauge couplings; 

• the different parameters can vary independently of each other. 



Conversely, the kinematic factors in the two models will be the same for a given set of masses 
and widths of the different Higgs bosons. Therefore, in those cases, many features of the signal, 

9 We use the same definitions of these couplings as in [5] . 
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Table 2: The maximal cross-sections in the 2HDM under consideration for Mh = 115 GeV, and 
Mn — 300, 500 and 700 GeV, respectively, with the following different cuts: inclusive, with 
primary cuts in eqs. ((HJ) — (EJ) , and with optimised cuts of eq. (jHJ) in the latter case also when 
only considering the H — ► hh resonant contribution. 

such as the differential distributions, will be similar to those of the MSSM even though the 
normalisation can be completely different. In fact, comparing Fig. H] with [5] through [7] we see 
that in the more general 2HDM the cross-sections can be more than two orders of magnitude 
larger than in the MSSM thus rendering a much larger potential for a detectable signal (as 
it will be discussed below). To be more quantitative on this we give in Tab. [2] the maximal 
inclusive cross-sections obtained in the scans for Mh = 115 GeV and Mh = 300, 500 and 
700 GeV. 

In order to study the potential signal in more detail we first of all apply the same primary 
cuts as in the case of the MSSM, those listed in eqs. (j5])-([7j). The resulting cross-sections 
are given Tab. El Comparing with the cross-section without the primary cuts we see that 
the reduction is substantial, but even so the signal cross-section can still be more than two 
orders of magnitude larger than in the MSSM scenario considered in subsection 13.1.11 and it is 
comparable to the background (see Tab. [T], specifically for low tan j3, where the Mh values in 
the two models are very similar). 

3.2 Signal-to-Background Differential Analysis 

In this section, we will continue the discussion of our numerical analyses limitedly to the Type 
II 2HDM considered so far. In order to enhance the statistical significance S/ \[B we studied 
several differential distributions for signals and background with the event selection of eqs. (jSJ)- 
([7]), with the aim of introducing optimised cuts, allowing at the same time to keep the signal 
event numbers at a reasonable level. To begin with, for simplicity, we have limited ourselves 
to use the contribution from the H — > hh resonance to the signal for Mh = 300 GeV in a 
scenario where the cross-section is close to maximal, with cos(/9 — a) = 1, Xnhh = 1000 GeV 
and Th = 30 GeV, when comparing with the background. 

The most sensitive distributions, able to discriminate between the signal and background, 
turn out to be the minimum transverse momentum of the forward/backward jets and the next- 
to-minimum invariant mass of the bb pairs, which we show in Fig. [81 (Although to a some more 



limited extent, also the minimum bb invariant mass is useful.) Before selecting a specific set 
up, we performed also a systematic analysis of the significance for different combination of cuts 
(30 GeV < mlf n < m^ ext - to - min < 100 GeV, 20 GeV < p^ d < 60 GeV). The best optimised 
cuts, on top of the basic ones of eqs. (|5])-([7j), that we found arH 

^fwd/bwd > 4Q GeV ^ ^ > 4Q GeV ^ m ncxt-to-min > g Q (g) 

We show in Fig. [9] the 46 invariant mass distribution for three signals (Mjj = 300, 500 and 
700 GeV with the widths Th = 30, 50 and 200 GeV, respectively) and the background after 
the optimised cuts of eq. (jHj) have also been imposed. For each of the three signals shown in 
the figure we have used the parameter space point which gives the maximal signal cross-section 
from the resonant H — > hh contribution when restricting the width to be less than 30, 50 
and 200 GeV, respectively. In this context we note that there are two effects which mainly 
determine the width of the signal distribution. On the one hand, the smearing of momenta we 
use gives a contribution to the measurable width of about 30 GeV. On the other hand, one of 
course has the intrinsic width of the H. 

Taking suitable mass windows around the peaks for the different Higgs mass values il- 
lustrated in Fig. El we obtain the maximal signal cross-sections, event numbers and statis- 
tical significances quoted in Tab. [H In order to calculate the signal cross-sections in the 
respective windows for different parameter space points, taking the actual width of the H 
into account, we rescaled the contribution from the H — > hh resonance with a factor cm h = 
(arctan [2(M# — mi)/Tj\ + arctan [2(jnu — Mji)/Yj\) where mi and mjj are the lower and up- 
per limits of the signal window, is the width of the signal distribution in parameter space point 



i estimated from Tj = J T 2 H _ + r| b with r 4b = 30 GeV being the width of the m 4b -distributions 
from finite detector resolution. (Notice then that cm h is a normalisation determined from sce- 
narios with = 30, 50 and 200 GeV for the different H masses.) Thus we approximate 
the cross-section in the window as cr pe ak = CM H &H^>hh- As a further requirement we also 
imposed that at least 50% of the signal cross-section after the optimal cuts comes from the 
H —>■ hh resonance such that the would-be-signal would not be obscured by other non-resonant 
contributions. 

The distributions of the signal cross-sections obtained in this way are given in Fig. [TUJ For 
illustration, the 5a limits at LHC, assuming an integrated luminosity of 300 fb _1 , o~ pca k > 2.7 
(2.3-3.3)0, 1.5 (1.3-1.8) and 0.8 (0.7-1.0) fb, for M H = 300, 500 and 700 GeV respectively, are 



10 Note that the efficiency of these is rather insensitive to the actual Higgs mass values, so that we have used 
the same set for any choice of the latter. 

11 The ranges given within parenthesis in this paragraph have been obtained by varying the factorisa- 
tion/rcnormalisation scale for the background by a factor of two around the default value, which makes the 
corresponding cross-section decrease by 30 % or increase by 50 % respectively. The reason for this is that, being 
essentially a six-jet cross-section, the background rate is proportional to a G s and it is therefore quite sensitive to 
the renormalisation scale. We also note that our default scale (p^ = Y^\Prjl n ) ^ as conservatively been chosen 
to be small so, if anything, our estimate of the final signal-to-background rates should be regarded as conserva- 
tive. In a real experiment one should of course attempt to use the sidebands for background normalisation. 
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Table 3: Number of events and significances for Mh = 115 GeV and Mh = 300, 500, 700 GeV in 
the respective best case scenarios, for a 46-tagging efficiency of (50%) 4 and after the optimised 
cuts of eq. (jHJ). The assumed integrated luminosity at LHC and SLHC are 300 fb -1 and 
3000 fb -1 , respectively. 

also illustrated and the fractions of parameter space points which gives cross-sections larger 
then this are 27 (24-30), 8 (5-12) and 0%, respectively. The corresponding numbers for the 
SLHC with 3000 fb" 1 are 43 (41-45), 31 (28-33) and 2 (0-2)%, respectively. Thus even at 
the SLHC we find no scope of observing a Mh = 700 GeV resonance in the channel under 
investigation. 

Finally we have also investigated the effects of restricting the allowed parameter space from 
tree-level unitarity by putting harder constraints on the maximal eigenvalue of the scattering 
matrices, A max . For this purpose, Fig. [TT] shows the signal cross-sections obtained in the scan as 
a function of A max . From the figure it is clear that the results (at least for Mh = 300 and 500 
GeV) are not sensitive to the precise value used for applying the unitarity constraint. On the 
other hand, applying a much harder constraint of the order A max ^4(12) (instead of A max < Sn) 
essentially leads to that the sensitivity for detection at the LHC is more or less washed out for 
Mh = 300 (500) GeV. The same also holds at the SLHC assuming an integrated luminosity of 
3000 fb" 1 . 

4 Conclusions 

We would like to conclude our paper by stating that, at both the LHC and SLHC, there exists 
a great potential to extract a H — > hh — > 46 resonance when Mh is constrained in the vicinity of 
115 GeV. This is a crucial result if one recalls that the detection of a sole Higgs resonance and 
consequent extraction of an Mh value may not point unambiguously to the underlying model 
of EWSB, not even in presence of further measurements of the heavier Higgs masses, M#, M& 
and/or M H ±. 

For example, the 2HDM considered here may be realised in a configuration wherein all 
visible Higgs masses are degenerate with those of the MSSM. Under these circumstances, we 
have proved that 

• it is not possible to extract an H — * hh — > 46 resonance from vector-boson fusion in the 
MSSM (not even if Mh is known) whilst 



• the opposite case is true in a substantial fraction of the parameter space of our 2HDM 
(even if Mh is not known), thereby enabling one to possibly measure the triple-Higgs 
coupling X Hhh . 

The latter is a Lagrangian term, which is different between these two models even when their 
patterns of Higgs masses and couplings to SM objects are the same, that would give a unique 
insight into the underlying EWSB mechanism. 

To be more specific our results show that in the most favourable scenario with Mh = 300 
GeV up to 27 (43) % of the parameter space would give a 5a signal at the (S)LHC assuming an 
integrated luminosity of 300 (3000) fb _1 when using the standard tree-level unitarity require- 
ment on the J = partial waves, Re(ao) < 1/2. These results are not sensitive to the precise 
value used for applying the unitarity constraint, albeit for very strong constraints the sensitivity 
for detecting the signal goes away. In the case of Mh = 500 GeV the fraction of parameter 
space probed is smaller with up to 8 (31) % giving a 5a signal, whereas for Mh = 700 GeV 
there is essentially no sensitivity at all. 

Despite we lack a full Monte Carlo simulation we believe to have incorporated the most 
critical aspects of the latter so that we do not expect more realistic studies (including parton 
shower, hadronisation, heavy hadron decays and detector effects) to affect too strongly our 
conclusions. 

Finally, we are currently pursuing other work along the directions outlined here, covering 
the case of lightest (neutral) Higgs boson pair production in the case of Higgs-strahlung and in 
association with heavy quarks [30J. 
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Diagrams by FEYNMANDRAW 




Figure 1: Feynman diagrams for q\q 2 — > qzqih^h^. Depending on the (anti)quark flavour combi- 
nation, the W ± - and Z-mediated graphs may not interfere. Besides, for final state (anti)quarks 
of different flavours, only half of the diagrams survive. 
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Figure 2: The masses of the neutral CP-even Higgs bosons as a function of the CP-odd one, 
for two choices of tan f3, assuming the MSSM in Maximal Mixing configuration. 
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Figure 3: The relevant couplings of the neutral CP-even Higgs bosons entering the production 
process in flTJ as a function of the CP-odd Higgs boson mass, for two choices of tan/5, assuming 
the MSSM in Maximal Mixing configuration. 
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Figure 4: The inclusive cross-sections (as denned in the text) for vector-boson fusion in (CQ), 
followed by hh — > bbbb decays, as a function of the CP-odd Higgs boson mass, for two choices 
of tan f3, assuming the MSSM in Maximal Mixing configuration. 



qq -> qqhh -> qqbbbb in 2HDM, m h =l 15 GeV, m H =300 GeV 
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Figure 5: The dependence of the inclusive cross-section qq^ — > qq^hh — > qq^bbbb in the 
2HDM under consideration on the different parameters when scanning over f 0000 parameter 
space points for M h = 115 GeV and M H = 300 GeV. (Note that M A and M H ± are free 
parameters.) 



qq -> qqhh -> qqbbbb in 2HDM, m h =l 15 GeV, m H =500 GeV 
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Figure 6: The dependence of the inclusive cross-section qq^ — > qq^hh — > qq^bbbb in the 
2HDM under consideration on the different parameters when scanning over f 0000 parameter 
space points for M h = 115 GeV and M H = 500 GeV. (Note that M A and M H ± are free 
parameters.) 



qq -> qqhh -> qqbbbb in 2HDM, m h =l 15 GeV, m H =700 GeV 
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Figure 7: The dependence of the inclusive cross-section qq^ — > qq^hh — > qq^bbbb in the 
2HDM under consideration on the different parameters when scanning over f 0000 parameter 
space points for M h = 115 GeV and M H = 700 GeV. (Note that M A and M H ± are free 
parameters.) 



> 

CD 
O 



6 

X) 

\ 

b 

T3 



2.5 



1.0 



I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 

solid: background n 
dash: signal 




75 100 125 

m bb (next— to— min) (GeV) 



> 

CD 

o 



Pi 



b 
T3 



' I ' 1 1 1 I ' 1 ' ' I 1 

solid: background 
dashed: signal 




I 



100 150 200 

p T (q)— min (GeV) 



Figure 8: The distribution of the next-to-minimum bb invariant mass (left) and of the minimum 
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tagging jet transverse momentum (right) for the signal (cross section for qq 
qq^bbbb in a close to best-case scenario for Mh = 300 GeV) and the background. The basic 
cuts of eqs. (jSI)-© are imposed. 



1 1 I ' 11 1 I ' 

solid: background 




300 400 600 800 1000 

m 4b (GeV) 

Figure 9: The differential cross-section da(qq^ —> qq^hh — ► qq^bbbb) / dm^ in the best case 
scenarios for Mh = 300, 500 and 700 GeV obtained when scanning over the available param- 
eter space restricting the width T H to be less than 30, 50 and 200 GeV, respectively. When 
calculating the signal distributions the actual widths have been assumed to be Th = 30, 50 
and 200 GeV, respectively. 
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Figure 10: Distributions of the resulting cross-sections gg^'-* — > gg^ ">hh -> qq n bbbb in the 2HDM 
under consideration using the optimal cuts obtained in a scan over 10000 parameter space points 
(the area is normalised to 1 for the cross-section with optimal cuts) for three different sets of 
Higgs boson masses as indicated in the respective plots. The solid line shows the results with 
optimal cuts, the dashed line shows the resonant contribution from the H — > hh processes and 
the long dashed line shows the resonant contribution in the respective signal windows requiring 
that at least 50% of the cross-section comes from the H —>■ hh resonance. The vertical line 
corresponds to the 5a-limit at LHC assuming 300 fb _1 and the integral of the curves to the 
right of it gives the percentage of parameter space points where the resonant cross-section is 
larger than this. 
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Figure 11: The distributions in resulting signal cross- sect ions gg^'- 1 — > qq^hh — > qq^bbbb in 
the 2HDM under consideration using the optimal cuts obtained in a scan over 10000 parameter 
space points as a function of the maximal eigenvalue A max of the scattering matrix for three 
different sets of Higgs boson masses as indicated in the respective plots. The upper (lower) 
horizontal line corresponds to the 5cr-limit at (S)LHC assuming 300 fb _1 (3000 fb _1 ). 



